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Nonresonant inelastic x-ray scattering (NIXS) is a possible tool to detect charge excitations in 
electron systems. In addition, multipole transitions at high-momentum-transfer region open a new 
possibility to determine orbital-selective electronic excitations in multi-orbital itinerant Sd elec¬ 
tron systems. As a theoretical example, we chose the antiferromagnetic state of iron arsenides and 
demonstrate that the orbital-selective excitations are detectable by choosing appropriate momentum 
transfer in NIXS. We propose that both NIXS and resonant inelastic x-ray scattering are comple¬ 
mentary to each other for fully understanding the nature of orbital excitations in multi-orbital 
itinerant electron systems. 
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I. INTRODUCTION 

When the incident photon energy of inelastic x-ray 
scattering is far from the absorption edge, the scattering 
occurs without resonance enhancement. This is called 
nonresonant x-ray scattering (NIXS). NIXS has been 
used to measure both phonon and electronic excitations. 
Among electronic excitations, charge excitations are ca¬ 
pable by NIXS, since photons are predominantly scat¬ 
tered by charge. Therefore, the dynamical density struc¬ 
ture factor is detectable with NIXS^i^ and thus NIXS 
has some complementary aspects with inelastic neutron 
scatterings^ 

The plane wave involved in the matrix element of NIXS 
can be expanded with respect to momentum transfer. 
In the high momentum-transfer region, multipole transi¬ 
tions become relevant, and thus NIXS using high-energy 
x-ray has a unique capability of measuring higher multi¬ 
pole transitions. By using this characteristic of NIXS, a 
possible observation of orbiton excitations in manganites 
has been proposeds^ Dipole-forbidden local d-d excita¬ 
tions in Mott insulators, NiO and CoO, have success¬ 
fully been observed via high momentum-transfer NIXS 
measurementsThe experimentally observed data have 
been analyzed theoretically on the basis of the first- 
principles calculations and cluster approachesFur¬ 
thermore, effective operators for the d-d excitations have 
also established and their application to manganites has 
been proposed^^ 

In itinerant 3 d electron systems, there is no clear local¬ 
ized d-d excitation but corresponding excitations related 
to orbital degrees of freedom exist. Therefore, it would 
be challenging to detect such orbital excitations by NIXS 
in high momentum-transfer region. A good example of 
multi-orbital systems where the orbital degree of free¬ 
dom is believed to be crucial for their physical properties 
would be iron-based superconductors 

In this paper, we demonstrate a possible realization 


of orbital-selective electronic excitations in the antifer¬ 
romagnetic phase of iron arsenides by NIXS. The an¬ 
tiferromagnetic phase is obtained by using the mean- 
field approximation for a five-orbital Hubbard model and 
dynamical susceptibilities are calculated within the ran¬ 
dom phase approximation (RPA). We find that specific 
orbital-selective excitations are detectable by choosing 
appropriate momentum transfers. Based on this finding, 
we propose that NIXS is a complementary experiment to 
resonant inelastic x-ray scattering (RIXS) tuned for Fe 
L absorption edge, from which the orbital characters of 
excitations can be extractedJ^ 

The rest of the paper is organized as follows. In Sec.HIl 
the model on the antiferromagnetic phase of iron-arsenide 
superconductors is introduced. In Sec. IIIII the theory of 
NIXS is briefly introduced. The matrix element of the 
plane wave with respect to Fe^+ iron is shown in Sec. IE 
together with the calculated NIXS spectra for BaFe2As2. 
Finally, a summary is given in Sec. [Vl 


II. MODEL HAMILTONIAN 

We start with a multi-band Hubbard Hamiltonian for 
irons in iron-based superconductors given by Hg = Hq J- 
Hj. The noninteracting Hamiltonian Hq is given by 

i,j cr,ij,,u 

where ct ^ ^ creates an electron at site i with orbital /i 
and spin a. Aij = Vi — Vj, where is the position 
of site i. and t{Aij] //, u) are the on-site energy and 
hopping integral, respectively. The on-site energies and 
the hopping integrals are taken from Kuroki et ah^ The 
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interaction Hamiltonian Hj is expressed asi^ 

i,fi 

-\- {U — 2 J) ^ 

f/-3J t t 

i,IJ,^u,a 

where U is the intra-orbital Coulomb interaction, and J 
is the Hand’s coupling. Here, we assume that the pair 
hopping is equal to J. 

In order to describe an ordered phase, we con¬ 
struct the mean-field Hamiltonian from 

and self-consistently solve mean-field equations con¬ 
taining the order parameter defined by = 

with the ordering vector Q, 
where N is the number of the lattice points, c]^ ^ = 

1 ^ and the average (• • •) is taken at 

zero temperature. For an antiferromagnetically ordered 
state, we can rewrite the sum of the wave vectors as 
^ Sfco Sm=o,i ^ ko-^mQ, where the sum of 
ko is over the magnetically reduced Brillouin zone (BZ). 

We consider a stripe-type antiferromagnetic order with 
Q = (7r,0) observed in the parent compound of iron- 
arsenide superconductor, BaFe 2 As 2 . We set U = 1.2 eV 
and J = 0.22 eV to yield a magnetic moment m = 0.8//^ 
{fiB is the Bohr magneton) at n = 6.0^ where n is the 
electron density and m = We 

note that m is chosen to be close to the measured value 
in BaFe 2 As 2 


III. NONRESONANT INELASTIC X-RAY 
SCATTERING 


where (l)^{r) is the Wannier wave function of orbital ji. 
Using and dynamical susceptibility 

the NIXS intensity reads 

I{q,oj) = + G)a\r{-qv - G) 

lJi,y A,r 

X (4) 


with 


k,k' 


uj — Hd + Eq -^r i6 


:X^k,iy,a^k-\-qr,fi,cr)‘ (5) 


In this paper, the dynamical susceptibilities are calcu¬ 
lated within RPA for the antiferromagnetically ordered 
phase given by 

The plan wave in ([3]) can be expanded with re¬ 
spect to the partial waves. 


oo I 

d<i- = Y, E ^T^i'3Mr)Yira{e„<t>,)Yil,{er,^r), (6) 

Z=0 m=—l 


where ji{x) is the spherical Bessel function, is 

the spherical harmonics, and and (j)q(r) ^im repre¬ 
sent polar angles of the q{r) vectors. When q is located 
within the reduced BZ, i.e., G = 0, the I = 0 component 
dominates the expansion (| 6 ]).— In this case, pq ([2]) can be 
approximated by the charge operator given by 

pQ ~ ^Qr — ^I,^,cr^^+gr,M,cr- (7) 

k,iJ,,a 

When q is larger than the reduced BZ, that is, G 7 ^ 0, 
non-zero I contributes to the expansion ©■ This means 
that NIXS has a unique capability of measuring higher 
multipole transitions. 


The spectral intensity of NIXS at zero temperature is 
given by 


IV. CALCULATED RESULTS AND 
DISCUSSIONS 


= - 


1 


u! — Hd + Eq + i5 


:xPq)^ ( 1 ) 


where Eq is the ground-state energy, (5 is a small positive 
value. By taking q = q^ ^ G with the reciprocal vector, 
G, and the reduced momentum in the reduced BZ, q^^ 
the operator pq is defined as 

Pq ~ + G) ^l,zy,cr^^ + 9r,Ai,Cr5 (2) 

k,a 

with the coefficient a^u{q) defined as 

anu{q) = J 4>l{r)E’^"^4>^,{r)dr, (3) 


In order to examine the q dependence of the coefficient 
([3|) for iron arsenides, we simply use hydrogen-like atomic 
Sd wave functions with the effective core charge Zeff = 8 
corresponding to Fe^+ ion, instead of the Wannier wave 
functions. We plot OL^y{q) along q/{27r) = [id, 0,0] 
in Fig. [Ha). In a parent compound of iron-arsenide 
superconductor BaFe 2 As 2 (Bal22), [1,0,0] corresponds 
to 27r/a = 27r/2.80 = 2.24A“^, where a is the Fe-Fe 
bond length which is the half of the low-temperature or¬ 
thorhombic lattice constant. When H is less than one, 
i.e., within the first BZ (G = 0), the diagonal components 
of a^i^{q) give almost the same value and there are very 
small off-diagonal components, leading to the charge op¬ 
erator ©• With increasing id, the diagonal components 
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FIG. 1. (Color online) The coefficient a^u{q) defined in (|3]) 
for Fe^^ ion. (a) q/{27r) = [if, 0,0]. (b) The orbital diagonal 
part of a^^{q) along the [0, 0, L]-[l/2, 0, L] direction for each 
value of L. 
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become unequal and off-diagonal one for // = 3z‘^ — 
and u = emerges as shown in Fig. [IJa). 

Figure [2][a) shows the calculated NIXS in the antifer¬ 
romagnetic phase of Bal22 along q = qr = 27r[0,0,0] to 
27r[l/2,0, 0]. We note that the calculation has been per¬ 
formed for a two-dimensional model given by Hd- In the 
figure, we use as the unit of energy to evaluate the 
excitation energy under the presence of renormalization 
effects. If there are no renormalization effects, Ej. = 1 eV, 
but we guess that the value of E^ is around 0.5 eV as dis¬ 
cussed separatelyThe spectra shown in Fig. [2][a) are 
almost equivalent to the dynamical charge susceptibility, 
and thus it is difficult to select orbital-resolved excita¬ 
tions. 

Since our purpose is to predict orbital-selective exci¬ 
tations by NIXS, it is crucial for identifying the region 
of q where only a specific orbital contribution remains in 
a^y{q). After examining a^y{q) in detail, we find such 
a q region where only two orbital diagonal components 
are dominant near G = 27r[0,0, 8]. To see this we show 
a^^{q) along the [0, 0,1/]-[l/2, 0, L] direction for a given 
integer L in Fig. [IJb). We note that in the interval 1 
to 2 in the horizontal axis, for example, the diagonal 
components along q^/i^Ti) = [0,0,0] to [1/2,0,0] with 
G/(27r) = [0,0,1] are plotted. Near L = 8, we can find 
that y = xy and y = x‘^ — y‘^ dominate a^^{q). This 
opens a possibility of observing orbital-selective excita¬ 
tions if one tunes G around this region. We note that 
along this momentum-transfer direction the orbital off- 
diagonal components are very small as compared with 
diagonal ones. 

As a demonstration, we plot in Fig. Hb) calculated 
NIXS for Bal22 along qr/{27r) = [0,0,0] to [1/2,0,0] 
starting from G/(27r) = [0,0,8]. We note again that the 
calculation has been done in the two-dimensional model 


FIG. 2. (Golor online) Galculated NIXS spectra for the anti¬ 
ferromagnetic phase of Bal22. (a) Along the [0, 0, 0]-[l/2, 0,0] 
direction, and (b) along the [0, 0, 8]-[l/2, 0, 8] direction. E^ is 
the energy unit introduced to take into account the renormal¬ 
ization of the band. 



FIG. 3. (Golor online) The coefficient a^^uiq) defined in (|3|) 
for Fe2+ ion. q/27v is along [H,H,S]-[H ^ 1/2, FT, 8]. Each 
panel corresponds to each orbital ly. 


and the z component of G is taken into account through 
the matrix element a^y{q). There appears a momentum- 
dependent strong intensity around u = O.SE'r, which is 
different from Fig. [TJa) where all of orbitals can con¬ 
tribute to the intensity. The intensity comes from either 
an xy to xy excitation or an x‘^ — y‘^ to x‘^ — y‘^ excita¬ 
tion. Analyzing the wave functions in the ground state 
and corresponding excited states, we find that the strong 
intensity comes from the xy-xy excitation. 

In multi-orbital systems such as Bal22, there are char¬ 
acteristic orbital off-diagonal excitations, some of which 
have been predicted by the previous study by the same 
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FIG. 4. (Color online) Calculated NIXS spectra for antifer¬ 
romagnetic phase of Bal22 along the [1,1, 8]-[3/2,1, 8] direc¬ 
tion. Er is the energy unit introduced to take into account 
the renormalization of the band. 
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FIG. 5. (Color online) The same as Fig. |4]but for the para¬ 
magnetic phase of Bal22. 


authors in connection with RIXS^^ In order to detect 
such orbital off-diagonal excitations, we again need to 
identify G where the off-diagonal components are en¬ 
hanced in We find that such a situation oc¬ 
curs along the [H^ 8] direction. Figure [3] exhibits 

a^y{q) along + 1/2, i7,8] for a given H. 

When u = 3z^ — and = 1, the magnitude of 
<^xy, 3 z‘^-r‘^{Q)\ is as large as the diagonal contribution 
Oi3z‘^-r‘^,3z‘^-r‘^{Q)\‘ Similarly, when u = x‘^ the off- 
diagonal components c^yz,x‘^-y‘^{Q) and oizx,x^-y^{Q) be¬ 
come large when H = 1. This implies possibility for the 
observation of orbital off-diagonal excitations in Bal22. 

In Fig. m the NIXS intensity from [1,1,8] to [3/2,1, 8] 
{H = Im Fig. [3]) is shown. Strong intensity at cj = 0.3F^r 
is clearly shown, whose origin is attributed to excitations 
from 3z‘^ — to xy by the simulated RIXSJ^ Figure [4| 
is, thus, a prediction to NIXS experiment in the antifer¬ 
romagnetic phase of Bal22 for detecting an orbital off- 
diagonal excitation from 3z‘^ — to xy. We note that 
the strong enhancement is a consequence of two contribu¬ 
tions: one is the sharp occupied and empty peaks with 
dominantly 3z‘^ — and xy components, respectively, 
in the density of statesi^ and the other is the enhanced 
\(^xy,3z‘^-r‘^{Q)\ shown in Fig. [3l 

It is interesting to compare the NIXS intensity at the 
same q with that in the paramagnetic phase of Bal22. In 



FIG. 6. (Color online) The difference intensity of NIXS be¬ 
tween the antiferromagnetic phase (Fig.[3|) and the paramag¬ 
netic phase (Fig. [5|). 

the paramagnetic phase, the calculated NIXS intensity is 
relatively broad as shown in Fig. [5] in contrast with the 
antiferromagnetic case in Fig. ID The broadness predom¬ 
inantly comes from less structured density of states in 
the paramagnetic state. As a result, the difference inten¬ 
sity shown in Fig. [6] enhances positive intensity around 
00 = 0.3Er, but it is negative below oo = 0.3Er. We thus 
propose to measure the spectral change of NIXS across 
the Neel temperature in Bal22. Although the scattering 
intensity calculated in the present study may be smaller 
than the phonon scattering by three or more orders of 
magnitude, as is the case of the d-d excitations in NiO^ 
we expect that the spectral features proposed would be 
observed in the near future like the d-d excitations in 
NiO. 


V. SUMMARY 

We have demonstrated theoretically orbital-selective 
electronic excitations in the antiferromagnetic phase 
of iron arsenides by employing NIXS. The antiferro¬ 
magnetic phase for Bal22 has been obtained by using 
the mean-field approximation for a five-orbital Hubbard 
model and dynamical susceptibilities have been calcu¬ 
lated within RPA. Using the formalism for NIXS, we have 
calculated its intensity for specific momentum-transfer 
regions where selected orbital contributions to NIXS are 
maximized. We have clearly shown that the specific 
orbital-selective excitations are detectable by choosing 
appropriate momentum transfers. For Bal22, we have 
predicted that the difference intensity along the specific 
momentum transfer picks up a specific orbital-selective 
excitation. Based on this finding, we propose that NIXS 
and RIXS are complementary to each other for identify¬ 
ing orbital excitations in itinerant 3d electron systems. 
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